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ABSTRACT 

Aims. We aim at characterizing the inward transition from convective to radiative energy transport at the base of the convective 
envelope of the solar-like oscillator HD 52265 recently observed by the CoRoT satellite. 

Methods. We investigated the origin of one specific feature found in the HD 52265 frequency spectrum. We modelled the star to 
derive the internal structure and the oscillation frequencies that best match the observations and used a seismic indicator sensitive to 
the properties of the base of the envelope convection zone. 

Results. The seismic indicators clearly reveal that to best represent the observed properties of HD 52265, models must include 
penetrative convection below the outer convective envelope. The penetrative distance is estimated to be ~ 0.95H F , which corresponds 
to an extent over a distance representing 6.0 per cents of the total stellar radius, significantly larger than what is found for the Sun. 
The inner boundary of the extra-mixing region is found at 0.800 ± 0.004 R where R = 1.3 R is the stellar radius. 
Conclusions. These results contribute to the tachocline characterization in stars other than the Sun. 
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1. Introduction 

Low-mass main-sequence (MS) stars have convective envelopes 
in which the chemical elements are mixed on short time scales 
and - at least in the deeper parts of the envelope - the energy 
transport by fluid elements can be treated as an adiabatic process. 
In the standard description, convective zones (CZ) are regions 
where the Schwarzschild criterion is fulfilled. Their boundaries 
lie at the border where the adiabatic temperature gradient equals 
the radiative one. It is expected, however, that fluid elements 
penetrate into the adjacent radiative zone due to their inertia (see 
e.g. Zahn 1991). In the Sun, the transition region (tachocline) is 
believed to be the site where the dynamo originates. 

Penetrative convection (PC) corresponds to efficient convec- 
tive heat transport - and material mixing - by downward flows 
that establish a close to adiabatic temperature stratification be- 
low which the downward plumes are no longer able to modify 
the temperature stratification that remains close to radiative. The 
extent of penetrative convection and/or overshoot in stars cannot 
be derived from first principles and is still largely unknown. 

A crucial point therefore is to find observational signatures 
of penetrative convection in low-mass MS stars. This can be 
achieved by means of asteroseismology. The abrupt change of 
energy transport from a convective to a radiative regime is visible 
in the sound speed profile and impacts the oscillations frequen- 
cies as well as some characteristic frequency spacings, which 
then show an oscillatory (periodic) behaviour (Gough 1990; 
Roxburgh & Vorontsov 1994). Owing to the importance of the 
tachocline region, helioseismic studies have attempted to mea- 
sure the extent of the PC in the Sun (Berthomieu et al. 1993; 
Monteiro et al. 1994). Christensen-Dalsgaard et al. (2011) have 



recently found that convective envelope overshoot is necessary 
over an estimated extent of Q31Hp (Hp is the pressure scale 
height). The high-quality solar data and wide available range of 
values of mode degrees t allowed Christensen-Dalsgaard et al. 
to also show that the transition between the convective and the 
radiative stratification must be smooth, intermediate between a 
classical - ballistic - overshoot formulation and a no-overshoot 
one. 

Solar-like oscillations have been identified in many stars first 
from the ground, then from space by the CoRoT (Baglin et al. 
2002) and Kepler (Koch et al. 2010) high-precision photometry 
missions, leading to an accuracy in frequency measurements of a 
few tenths //Hz. Theoretical studies of the effects on the oscilla- 
tions frequencies of PC at the base of a CZ have been conducted 
using low-degree modes, the only ones expected to be detectable 
in stars (see e.g. Roxburgh 2009, and references therein). The 
period of the oscillatory component is found to be related to the 
location of the discontinuity inside the star and its amplitude to 
depend on the height of the discontinuity. One then expects that 
the amplitude of the oscillatory signal grows with an increasing 
extent of PC that causes a larger jump from the adiabatic temper- 
ature gradient to the radiative one below. Moreover, for a larger 
PC extent, the discontinuity is located deeper inside the star and 
the period of the oscillation is expected to be longer (Roxburgh 
1993). 

Ballot et al. (201 1) analysed the CoRoT oscillation spectrum 
of the star HD 52265 (HIP 33719), a high-metallicity GOV star 
hosting an exoplanet. They found a typical p-mode solar-like 
spectrum, and identified 28 reliable low-degree p-modes of de- 
grees £ — 0, 1 , 2 and order n in the range 14-24 (see their Table 
4). The frequencies v„/ are in the range 1500 - 2550 //Hz with a 
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frequency at maximum amplitude v max = 2090 + 20yuHz. The er- 
ror on each frequency is a few tenths of yt/Hz. Such a high quality 
data set enables to probe the interior structure of the star. Here 
we report on one evidence for penetrative convection below the 
upper convective region of HD 52265 as indicated by its internal 
structure modelling. 

2. HD 52265: observations and modelling 

2.1. Global and seismic observational constraints 

To model HD 52265 we adopted the effective temperature 
r eff = 6120 + 110 K and metallicity [Fe/H] = 0.22 + 0.05 
dex that we derived from an average of 20 spectroscopic de- 
terminations reported since 2001. We adopted the luminosity 
L = 2.053 + 0.053L o that we derived from the Hipparcos paral- 
lax (van Leeuwen 2007), Tycho magnitude and bolometric cor- 
rection calculated according to VandenBerg & Clem (2003). We 
related the ratio of heavy elements mass fraction Z to hydro- 
gen mass fraction X to [Fe/H] through [Fe/H] = log(Z/X) - 
log(Z/X) Q and adopt (Z/X) Q = 0.0244 from the Grevesse & 
Noels (1993) solar mixture. 

In the following we will focus on the ability of models to 
reproduce the values of the seismic indicators rroi(n) and rrio(n) 
introduced by Roxburgh & Vorontsov (2003), which are defined 
as 

rr i(ri) = ddoi(ri)/Avi(ri) ; rr w (n) - dd w (n)/Av (n), (1) 
where 

ddoi(n) = -(v„_i j0 - 4v„_i j i + 6v„ j0 - 4y„,i + y„+i j0 ) (2) 

o 

ddi (n) = -^(v„-i,i - 4v„, + 6v n>/ - 4v„ +1 , + v n+u ), (3) 

where Av((n) = v„+i,f - v„/ are the standard large frequency 
separations (Tassoul 1980). Roxburgh & Vorontsov showed that 
the ratios rroi and rrm (hereafter rroi/io) are sensitive to the 
sharp variation of the sound speed in the transition region be- 
tween the upper CZ and the radiative layers beneath (see also 
Otf Floranes et al. 2005). Roxburgh (2009), using high-quality 
solar data, confirmed that these indicators exhibit an oscilla- 
tory behaviour as a function of the frequency. From the 28 re- 
liable individual frequencies determined by Ballot et al. (201 1), 
we calculated the fToi/io(«) ratios. The variations of rroi(n) 
(resp. rrio(n)) as a function of frequency v„ ( o (resp. v„ i) are 
plotted in Fig 1. Ballot et al. already noticed a periodic varia- 
tion of the separations do\(n) = Vo,«-i - 0.5(vi „ + vi „_i) and 
dw(n) = -Vi, n + 0.5(vo in + vo,n+i) and suggested that it could be 
the signature of the CZ boundary. The periodic variation is even 
clearer in the rroi/io(«) separation, and Section 3 is devoted to 
confirming the origin of this periodic signal and to characteriz- 
ing it. 

2.2. Internal structure models and their oscillations 

We modelled HD 52265 with the evolution code Cesam2k 
(Morel & Lebreton 2008) and the input physics and param- 
eters described in Lebreton (2011). We considered the model 
for penetrative convection below the convective envelope pro- 
posed by Zahn (1991). In this model the distance of fluid pen- 
etration into the radiative zone reads L p = (£/xp)Hp, where 
Xp = (31og^/51ogP) a d is the adiabatic derivative with respect 
to pressure P of the radiative conductivity^ = 16crr 3 /(3p/c) 
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Fig. 1. Observed ratios rroi/io(«) (black squares) as a function of fre- 
quency v„ Q/ „ i for HD 52265 compared to the rr i/io(n) ratios in the 
models without PC = 0.0, magenta) and with PC = 1.3, red). 



(T,p,K,cr are the temperature, density, opacity, and Boltzmann 
constant, respectively). The free parameter £ is of the order 
of unity but has to be calibrated by comparing stellar models 
to observations. The seismic indicators rroi/io that we consid- 
ered here are sensitive to the change in the temperature deriva- 
tive hence to the transition from an unstable to a stable strat- 
ification. The amplitude of the periodic signal is smaller for a 
smoother transition and cannot be detected if the transition is 
highly smooth. We therefore concentrated on determining the 
adiabatic extent of the overshoot region. We imposed accord- 
ingly that the temperature gradient in the overshooting zone is 
the adiabatic gradient. 

Equilibrium models were calculated and adjusted to sat- 
isfy the constraints provided by the global parameters (L, T e s 
and surface [Fe/H]) and the 28 reliable observed frequencies of 
Ballot et al. (201 1). The model frequencies were calculated with 
the LOSC adiabatic oscillation code (Scuflaire et al. 2008) for 
the whole range of observed mode orders and degrees, and the 
observed and modelled seismic indicators were derived consis- 
tently. We corrected the model frequencies from the so-called 
near surface effects using the empirical correction proposed by 
Kjeldsen et al. (2008). The frequency differences may be sen- 
sitive to surface effects, however, we recall that Roxburgh & 
Vorontsov (2003) demonstrated that the ratios of small to large 
separations are quite independent of the surface treatment and 
are therefore efficient probes of the interior. 

We used the Levenberg-Marquardt minimization method as 
described in Miglio & Montalban (2005) to adjust the age and 
mass, the initial helium and metallicity, the mixing-length pa- 
rameter for convection, the convective core overshooting extent, 
and surface effects parameter so that the model of HD 52265 
fits the observations best, within the error bars. The goodness 
of fit was evaluated through minimization of the reduced x 2 '- 

X 2 = W>bs - IT 1 • 2^T (Kmod - ■Xi,ob s )/cr i , obs ) 2 , where 7V obs is 
the number of observational constraints considered, JCi, mo d and 
Xi bs are the modelled and observed values of the z™ parameter, 
respectively, and cr i obs is the error on Xi j0 b s . We point out that 
the convective core overshooting extent derived from prelimi- 
nary model calibrations of HD 52265 is low, lower than 0. \QH p . 
Because the properties of the central region do not affect the ef- 
fect we study at the transition region, the models presented here 
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Fig. 2. Top: Temperature gradient as a function of radius in a model with 
penetrative convection (continuous line) and a standard model (dashed 
line) of HD 52265. Bottom: Sound speed gradient with respect to the 

acoustic depth r(r) = J dr/c s for the same models. 



do not include core overshooting. The detailed properties of the 
models will be presented in a forthcoming paper. 



3. Extent of penetrative convection 

The best-fit models yield totals-values in the range 2-3 what- 
ever the value of the PC and therefore cannot be distinguished 
on that basis. The models have an age A ^2 Gyr, a mass 
M =s 1 .25 M Q and a radius R 1 .3 R Q in quite reasonable agree- 
ment with the values derived by Escobar et al. (2012) in their 
recent modelling of the star. We now examine the impact of PC 
on the rroi/io(«) frequency ratios. 

We first compared models with £ = 0.0 and 1.3. The varia- 
tions of the ratios ttoi/io(«) with frequency are shown in Fig.l. 
The amplitude of the oscillatory component is significantly 
smaller for £ = 0.0 than for £ = 1.3. The impact of a nonva- 
nishing £ on the structure of the model with £ = 1.3 compared 
with a model without PC is illustrated in Fig. 2. Differences ap- 
pear in the transition region below the convection zone while ev- 
erywhere else, the structures coincide. The top of Fig. 2 shows 
the behaviour of the temperature gradient below the CZ in the 
model of HD 52265. Without PC the CZ is located at r Sc , the 
Schwarzschild radius. The temperature gradient drops from its 
adiabatic value 0.4 at the border of the convective envelope and 
takes the value of the radiative gradient. In contrast, when PC 
takes place below rs c as in model with £ = 1.3, a discontinuity 
in the temperature gradient occurs at the bottom of the mixed 
region (including the PC extent), resulting in a nearby density 
discontinuity. This sharpens the variation of the adiabatic sound 
speed c s = (FiP/p) 1 ^ 2 at the base of the convective envelope in 
models including PC (Fj is the first adiabatic index). The bottom 
of Fig. 2 shows the sound speed gradient as a function of the 

rR 

acoustic depth r(r) = J r dr/c s . Rapid variations in the sound 
speed gradient can be seen, which correspond to the Hel ioniza- 
tion (r ~ 620 s) and to the inner limit of the external CZ. At the 
base of the CZ, the variation of dcjdr for the model including 
PC behaves as a near discontinuity (t ~ 2780 s), whereas is re- 
mains continuous in the standard model where only d 2 c s /dT 2 is 
marginally discontinuous (r ~ 2450 s). 



HD52265 
| = 1.3 
sinusoidal fit 
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Fig. 3. Difference of the ratios rroi/io(rc) as a function of frequency be- 
tween the standard model (no PC, £ = 0.0) and (i) the model with PC 
with an extent of £ = 1.3 (red continuous line), (ii) the observations 
(black). The dashed red line is a sinusoidal fit of the red curve (see 
text). Note that to calculate the differences, models were splined on the 
observed frequency grid. A detailed comparison with observations is 
provided in Fig. 4. 

Table 1. Fit of the signal S = Arr i /W (n) = rr 01/w (n) s - rr Q{no (n\ =0 . 
Case 1: Arr oi/lo (n) = (a/v + b/v 2 )cos(47tvT + 2</))+offset. Case 2: 
A/T i/io(n) = (a/v)cos(47rvr + 20) + (b/v 2 ) sin(47rvr + 2i/r)+offset. 
Model fit (M( =l 3 ) ; fit of raw observations (O) ; spline adjustment on 
observations (Os). 
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To examine the ability of models to reproduce the observed 
rroi/io oscillatory trend, the difference between the signals from 
the observations and models was measured with a specific x 1 
quantity where the variables are x- t = rroi/io(«). Considering 
models with increasing we found that the variation of x 2 
shows a minimum between £ = 1.2 and 1.3. The signals from 
models with £<l.lor£>1.4 lead to significantly larger^ 2 , 
excluding these models as acceptable. We therefore found a PC 
extent of f 1.25 + 0.10 for HD 52265. This translates into an 
overshoot distance of d m = 0.95 + 0.08 H p corresponding to 
d m = 0.060 + 0.004 R. 



4. Seismic location of the base of the adiabatic 
temperature stratification 

The structure of inner regions below the transition layer is un- 
affected by a change in It is therefore reasonable to assume 
that the long-term trend seen in Fig.l for both signals with and 
without PC is the same and has its origin deeper in the star. This 
enables us to remove it to isolate the oscillatory component aris- 
ing from the base of the CZ of model with £ = 1.3. We interpo- 
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lated the rroi/io(«) signals of the models with £ = 0.0 and 1.3 
on the same frequency grid and computed their differences. The 
differences plotted in Fig. 3 are clearly generated by including 
PC in the model with £ = 1.3. We obtain an almost pure sinu- 
soid, the characteristics of which can be estimated and related to 
the properties of the region of penetrative convection. Slightly 
different expressions of the sinusoidal component arising from 
the CZ base that is expected to be visible in several seismic in- 
dicators of distant stars have been established theoretically (see 
e.g. Roxburgh & Vorontsov 1994; Monteiro et al. 2000; Basu 
et al. 2004; Mazumdar 2005; Verner et al. 2006; Houdek & 
Gough 2007; Roxburgh 2009). Guided by them, we assumed that 
the model oscillation in Fig. 3, which stems from the difference 
Arr i/io(«) = ""oi/io(«)f=i.3 - ""oi/io(«)f=o.o, can take two nearly 
equivalent forms expressed in the caption of Table 1 . Roxburgh 
(2009) expressed the periodicity of the signal as 1 /(2T^ = i.3), 
which provides the acoustic radius T^r = i 3 = T,-r(r) at the base of 
the CZ. r(r) is the acoustic depth and r, = r{r = 0) = l/(2(Av)), 
where (Av) is the mean asymptotic large frequency separation. 
The amplitude is a slowly varying function of v with a term in 
v (resp. v~ 2 ), arising from a discontinuity in the first (resp. sec- 
ond) derivative of the sound speed. We list in Table 1 the values 
of Tf=i.3 and other parameters of the fit of the sinusoidal signal 
of Fig. 3 (xL =1.6 10~ 4 ). The results weakly depend on the as- 
sumed form for Arroi/io(«). We find Tf=i.3 - 2167 + 15 s, which 
is higher than values obtained by Monteiro et al. (2000) because 
Monteiro et al. considered less overshoot and slightly smaller 
stellar masses. These T^ = ; 3 values can now be compared with 
the theoretical value, T ov = r t — r(r ov ), directly computed from 
the equilibrium model with £ = 1.3, where r ov is the radius at the 
base of the adiabatic region and t, — 4913 s for the model with 
£ = 1.3. The fitted value T f=1 . 3 falls within ~ 35 s (~ 2%) of the 
resulting theoretical value T ov = 2132 s, the shift T^=i .3 - T ov 
being similar to that found for the Sun by Roxburgh (2009). 

The same procedure was applied to the observed signal, 
which is compared to the model with ^ = 0.0 and the sinusoidal 
fits are plotted in Fig. 4. The results of the fits (yjL = 10~ 3 ) are 
listed in Table 1. The period r o b s = 2100 ± 100 s agrees quite 
well with that derived for model £ = 1.3. We conclude that the 
base of the adiabatically stratified region of HD 52265 is located 
within 2% of the radius r m /R = 0.800 ± 0.004. 

5. Conclusions 

We have detected an oscillatory signal with a significant am- 
plitude in the variation of the separations rroi/io(«) with fre- 
quency for the C0R0T star HD 52265. A comparison with the 
same signal arising from appropriate stellar models shows that it 
cannot be reproduced unless penetrative convection is included 
at the base of the outer convective envelope. This is the first 
time that such a feature is firmly detected in a star other than 
the Sun. A best fit of the signal provides a measure of the ex- 
tent of the mixed region below the CZ, in terms of the proxy 
d m = 0.95 + 0.08 Hp. 

The periodic signal for HD 52265 is more pronounced than 
for the Sun (Christensen-Dalsgaard et al. 2011), with a longer 
period, and so is the measured extent of PC in terms of H p and 
normalized radius. We point out that HD 52265 is similar to the 
Sun in all aspects except for the higher metallicity. Therefore, to 
understand the amplitude difference between the two stars, it is 
important to investigate the impact of metallicity on the structure 
and dynamics of the tachocline. Progress is expected in a near 
future when seismic missions will provide similar observations 
for many stars spanning a wide metallicity range. 
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Fig. 4. Difference of the ratios rroi/io(i) as a function of frequency be- 
tween the standard model (no PC, £ = 0.0) and the raw observations 
(black) or a spline fit of observations (continuous blue line). The dashed 
lines are the fits of the raw and splined data (see Table 1). 



With the above results, we provided an additional clue that 
the physical description of the turbulent convection must be im- 
proved in stars, particularly at interfaces with radiative regions. 
Future modelling of the dynamics in the region of the tachocline 
will have to comply with our findings as well as with the results 
previously obtained for the Sun. 
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